%I made a change.

% Checks that the Fp and Fc values are reasonable

theta = linspace(0, pi, 20);
phi = linspace(-pi, pi, 20);
psi = linspace(0, pi, 20);

detNames = {'LHO', 'LLO', 'Virgo', 'INDIGO'};
for k=1:4;
    Det = LoadDetectorData(detNames{k});
    [Fp Fc] = ComputeAntennaResponse(phi, theta, psi, Det.d);
    max(Fp)
    max(Fc)
%     figure(k);
%     plot(psi, Fc)
%     figure(k+4);
%     plot(psi, Fp)
end    
